Broadening of a nonequilibrium phase transition by extended structural defects 
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We study the effects of quenched extended impurities on nonequilibrium phase transitions in the 
directed percolation universality class. We show that these impurities have a dramatic effect: they 
completely destroy the sharp phase transition by smearing. This is caused by rare strongly coupled 
spatial regions which can undergo the phase transition independently from the bulk system. We 
use extremal statistics to determine the stationary state as well as the dynamics in the tail of the 
smeared transition, and we illustrate the results by computer simulations. 
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In nature, thermal equilibrium is more of an excep- 
tion than the rule. In recent years, phase transitions 
between different nonequilibrium states have become a 
topic of great interest. A prominent class of nonequilib- 
rium phase transitions separates active, fluctuating states 
from inactive, absorbing states where fluctuations cease 
entirely. These absorbing state transitions have applica- 
tions ranging from physics to chemistry and to biology 
B lj| The generic universality class for absorbing 
state transitions is directed percolation (DP) 5]. Accord- 
ing to a conjecture by Janssen and Grassberger 0, all 
absorbing state transitions with a scalar order parameter, 
short-range interactions and no extra symmetries or con- 
servation laws belong to this class. Examples include the 
transitions in the contact process Q, catalytic reactions 
H, interface growth [jj, or turbulence However, 
despite its ubiquity in theory and simulations, clearcut 
experimental realizations of the DP universality class are 
strangely lacking . The only verification so far seems 
to be found in the spatio-temporal intermittency in fer- 
rofluidic spikes [T^ . 

A possible reason for this apparent discrepancy are im- 
purities, i.e., quenched spatial disorder. According to 
the Harris criterion \13L Il4| , the DP universality class is 
unstable against disorder, because the (spatial) correla- 
tion length exponent v± violates the inequality v±_ > 2/d 
for all spatial dimensionalities d < 4. Indeed, in the 
corresponding field theory, spatial disorder leads to run- 
away flow of the renormalization group (RG) equations 
Ra, destro ying the DP behavior. Several other studies 
|l6l IrH fl8l Il9| agreed on the instability of DP against 
spatial disorder, but a consistent picture has been slow to 
evolve. Recently, Hooyberghs et al. applied the Hamil- 
tonian formalism |2fJ to the contact process with spa- 
tial disorder 1211 . Using the Ma-Dasgupta-Hu strong- 
disorder RG 123 these authors showed that the transition 
(at least for sufficiently strong disorder) is controlled by 
an exotic infinite-randomness fixed point with activated 
rather than the usual power-law scaling. In many real 
systems, the disorder does not consist of point defects 
but of dislocations, disordered layers, or grain bound- 
aries. The effects of such extended defects are generically 
stronger than that of uncorrelated disorder, as has been 
shown by detailed studies of equilibrium systems rang- 



ing from the exactly solved McCoy- Wu model [2I] and 
several RG studies [H IM IH Hi to the discovery of 
infinite-randomness critical behavior in the correspond- 
ing quantum Ising model [23 ■ 

In this paper, we show that extended defects have an 
even more dramatic effect on non-equilibrium phase tran- 
sitions in the DP universality class; they destroy the 
sharp transition by smearing. This is caused by phe- 
nomena similar to but stronger than the usual Griffiths 
effects |3 IIH : rare strongly coupled spatial regions can 
undergo the transition independently from the bulk sys- 
tem. In the tail of the smeared transition, the spatial den- 
sity distribution is very inhomogeneous, with the average 
stationary density and the survival probability depending 
exponentially on the control parameter. The approach of 
the average density to this exponentially small stationary 
value occurs in two stages, a stretched exponential decay 
at intermediate times, followed by power-law behavior at 
late times. In the following, we derive these results for a 
disordered contact process, illustrate them by computer 
simulations, and discuss their generality and importance. 

Our starting point is the clean contact process 0, a 
prototypical system in the DP universality class. It is 
defined on a rf-dimensional hypercubic lattice. Each site 
r can be vacant or active, i.e, occupied by a particle. 
During the time evolution, particles are created at vacant 
sites at a rate An/ (2c?) where n is the number of active 
nearest neighbor sites and the 'birth rate' A is the control 
parameter. Particles are annihilated at unit rate. For 
small A, annihilation dominates, and the absorbing state 
without any particles is the only steady state. For large A 
there is a steady state with finite particle density (active 
phase). Both phases are separated by a nonequilibrium 
phase transition in the DP universality class at A = A°. 

We introduce quenched spatial disorder by making the 
birth rate A a random function of the lattice site. Ex- 
tended impurities can be described by disorder perfectly 
correlated in d c dimensions, but uncorrelated in the re- 
maining d r = d — d c dimensions. A is thus a function of 
r r which is the projection of the position vector r on the 
uncorrelated directions. For definiteness, we assume that 
the A(ry) have a binary probability distribution 

P[A(r r )] = (l-p) S[X(r r ) - A] + p S[X(r r ) - cA] (1) 
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where p and c are constants between and 1. In other 
words, there are extended impurities of density p where 
the birth rate A is reduced by a factor c. 

Let us now consider the effects of rare disorder fluc- 
tuations in this system. In analogy to the Griffiths phe- 
nomena [tII^ . there is a small but finite probability for 
finding large spatial regions devoid of impurities. These 
rare regions can be locally in the active phase, even if 
the bulk system is still in the inactive phase. For the 
largest rare regions this starts to happen when A crosses 
the clean critical point A^. Since the impurities in our 
system are extended, each rare region is infinite in d c di- 
mensions but finite in the remaining d r dimensions. This 
is a crucial difference to systems with uncorrelated disor- 
der, where the rare regions are finite. In our system, each 
rare region can therefore undergo a real phase transition 
independently of the rest of the system. Thus, those rare 
regions that are locally in the ordered phase will have a 
true nonzero stationary density, even if the bulk system 
is still in the inactive phase. 

The resulting global phase transition is very different 
from a conventional continuous phase transition, where 
a nonzero order parameter develops as a collective effect 
of the entire system, accompanied by a diverging corre- 
lation length in all directions. In contrast, in our system, 
the order parameter develops very inhomogeneously in 
space with different parts of the system (i.e., different 
r r regions) ordering independently at different A. Cor- 
respondingly, the correlation length in the uncorrelated 
directions remains finite across the transition. This de- 
fines a smeared transition. Thus, extended impurities 
lead to a smearing of the DP phase transition. 

We now use extremal statistics to derive the properties 
in the tail of the smeared transition, i.e., in the parame- 
ter region where a few active rare regions exist, but their 
density is so small that they can be treated as indepen- 
dent. We start with the stationary behavior. The prob- 
ability w for finding a rare region of linear size L r devoid 
of impurities is, up to pre-exponential factors, given by 



exp(-pi^) 



(2) 



with p = — ln(l — p). As discussed above, such a region 
undergoes a true phase transition to the active phase at 
some \ c {L r ) > \®. According to finite-size scaling [30l |. 



Ac(-Ct) — A c 



all 



(3) 



where <f> is the clean (d-dimensional) finite-size scaling 
shift exponent and A is the amplitude for the crossover 
from a d-dimensional bulk system to a 'slab' infinite in d c 
dimensions but but finite in d r dimensions. If the total 
dimensionality d = d c + d r < 4, hyperscaling is valid and 
<p = l/v± which we assume from now on. Combining 
@ and we obtain the probability for finding a rare 
region which becomes active at A c as 



(4) 



for A c - A° -> 0+. Here B = pA drU± - . The total (average) 
density p at a certain A is obtained by summing over all 



active rare regions, i.e., all regions with A c < A. Since the 
functional dependence on A of the density on any given 
active island is of power-law type it does not enter the 
leading exponentials but only the pre-exponential factors. 
Thus, the stationary density develops an exponential tail, 



p(X)~exp(-B(\-\° c )- d ^) , 



(5) 



reaching to the clean critical point A^. Analogous argu- 
ments can be made for the survival probability P(X) of 
a single seed site. If the seed site is on an active rare 
region it will survive with a probability that depends on 
A with a power law. If is not on an active rare region, 
the seed will die. To exponential accuracy the survival 
probability is thus also given by JSJl. 

The local (spatial) density distribution in the tail of 
the smeared transition is very inhomogeneous. On active 
rare regions, the density is of the same order of magni- 
tude as in the clean system. Away from these regions it 
decays exponentially. The typical local density pt yp can 
be estimated from the typical distance of any point from 
the nearest active rare region. From (JIJ we obtain 



exp [B(\ - A°)" 



/d r 



(6) 



At this distance, the local density has decayed to 
Ptyp ~ e- rt -/« ~ exp {-Cexp [B(X - X° c )- d ^/d r ] } 

where £o is the bulk correlation length (which is finite 
and changes slowly across the smeared transition) and 
C is a constant. A comparison with (JSJ shows that the 
relation between the typical and the average density is 
exponential, |logptyp| ~ p~ 1 ' dr , indicating an extremely 
broad local density distribution. 

We now turn to the dynamics in the tail of the smeared 
transition. The long-time decay of the density is domi- 
nated by the rare regions while the bulk contribution de- 
cays exponentially. According to finite size scaling |30| . 
the behavior of the correlation time £t of a single rare 
region of size L r in the vicinity of the clean bulk critical 
point can be modelled by 



--!/«/ j. 



(8) 



Here A = A — A^ > 0, z is the d-dimensional bulk dynam- 
ical critical exponent, and and z are the correlation 
length and dynamical exponents of a d r -dimensional sys- 
tem. Let us first consider the time evolution of the den- 
sity at A = A°. For A = 0, the correlation time © 
simplifies to £t ~ To exponential accuracy, the time 
dependence of the average density is obtained by 



P(t) 



J dL r exp (-pLf - Dt/Ll) 



(9) 



where D is a constant. Using the saddle point method to 
evaluate this integral, we find the leading long-time decay 
of the density to be given by a stretched exponential, 

_j.d r /{d r +z) 



lnp(t) 



(10) 
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FIG. 1: Left: Overview of the steady state density of a clean (p — 0) and a diluted (p — 0.3) system. Center: Logarithm of the 
density as a function of (A — A^) -1 / 2 for several dilutions p and L = 10 4 . The data are averages over 100 disorder realizations. 
The solid lines are fits to 10 with d r v±_ — 1/2. Right: Decay constant B as a function of — ln(l — p). 



For A > A°, we repeat the saddle point analysis with 
the full expression JSJ for the correlation length. For in- 
termediate times t < t x ~ (A — \^-{d r +z)v± ^ e Qi eca y f 
the average density is still given by the stretched expo- 
nential (|lf)|> . For times larger than the crossover time t x 
the system realizes that some of the rare regions are in the 
active phase and contribute to a finite steady state den- 
sity. The approach of the average density to this steady 
state value is characterized by a power-law. 

(11) 

The value of tp cannot be found by our methods since it 
depends on the neglected pre-exponential factors. 

We now illustrate the smearing of the phase transition 
by results of a computer simulation of a 2d contact pro- 
cess with linear defects (d c = d r = 1). To reach the 
rather large system sizes necessary to observe exponen- 
tially rare events, we consider a version of the contact 
process with infinite-range couplings in the correlated di- 
rection (parallel to the impurities) but nearest-neighbor 
couplings in the uncorrelated direction (perpendicular to 
the defect lines). While this infinite-range model will not 
be quantitatively comparable to a short-range contact 
process, it provides a simple example for the smearing 
mechanism. Moreover, since the smearing only relies on 
a single rare region undergoing a true phase transition, 
we expect that the results will qualitatively valid for a 
short-range contact process, too (with the appropriate 
changes to the exponents). 

Because the couplings in the correlated direction are 
of infinite range, this dimension can be treated exactly 
in mean-field theory. This leads to a set of coupled local 
mean- field equations for the local densities p x , 

d X(x) 

— p x = -p x + — — (1-Px) {Px-i + 2px + Px+i) ■ (12) 

These equations can be easily solved numerically. We 
study systems with several dilutions p = 0.2 . . .0.6 and 
sizes of up to L — 10 6 in the uncorrelated direction; the 
impurity strength is c = 0.2 for all calculations. 



To determine the stationary state we solve the equa- 
tions -§ip x — in a self-consistency cycle. Our results are 
summarized in Fig.^ The left panel shows the total (av- 
erage) density p for a clean (p = 0) and a dirty (p = 0.3) 
system. The clean system has the expected sharp phase 
transition at A = Aj? = 1 with the mean-field critical ex- 
ponent (3 — 1. The data of the disordered system seem 
to suggest a transition at A k, 1.04. However, a closer in- 
spection shows that the singularity is smeared. Note that 
the density data are essentially size-independent. There- 
fore, the observed rounding cannot be due to finite-size 
effects. We conclude that the smearing is an intrinsic 
effect of the infinite system. For comparison with the an- 
alytical prediction (JSJ, the central panel shows the loga- 
rithm of the total density as a function of (A — A°) -1 / 2 
for different impurity concentrations p. Note that in our 
infinite-range model u± = 1/0 = 1/2. The data follow 
eq. (jSJ over several orders of magnitude in p. Fits of the 
data to (J5J) are used to determine the decay constants B. 
The right panel of Fig^shows that these decay constants 
depend linearly on p = — ln(l — p), as predicted. 

To study the time evolution we numerically integrate 
the local mean-field equations (|12fl . starting from a ho- 
mogeneous initial state with p = 1. Fig. [21 summarizes 
our results for a system of size L = 10 6 with dilution 
p = 0.5. The main panel shows a log-log plot of lnp 
vs. t. This allows us to test the stretched exponential 
behavior predicted in ljTU|) for the time dependence of p 
at the clean critical coupling A = A° = 1. We find that 
the data indeed follow a stretched exponential with an 
exponent of approximately 0.32 in excellent agreement 
with the analytical prediction d r /(d r + z) = 1/3. For 
A > A°, the decay initially follows the stretched expo- 
nential, but eventually the density approaches its finite 
steady state value. The inset of Fig. [3 shows a log-log 
plot of p(t) — Pstcady vs. t. The data clearly display 
power-law behavior in agreement with (|11|) . A fit to this 
equation gives an exponent value of tp « 2.6. 

To summarize, we have demonstrated that extended 
impurities destroy the sharp DP phase transition in the 
contact process by smearing. In the remaining para- 
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FIG. 2: Density p vs. time t for a system of size L = 10 , 
dilution p — 0.5 and several A (averages over 25 disorder re- 
alizations). Solid line: Fit of the A = 1 data (t > 100) to eq. 
1101 1 giving an exponent of approx. 0.32. Inset: Approach to 
the steady state density for A = 1.01. Solid line: Fit of the 
data for t > 100 to 1111) . giving an exponent of tp ~ 2.6. 

graphs we discuss the generality of our findings as well as 
their relation to Griffiths effects [14L |29| . The origins of 
Griffiths effects and the smearing found here are very sim- 
ilar, both are caused by rare large spatial regions which 
are locally in the ordered phase. The difference between 
them is a result of disorder correlations. For uncorrelated 
disorder, the rare regions are of finite size and cannot 
undergo a true phase transition. Instead, they fluctuate 
slowly. In contrast, if the rare regions are infinite in at 
least one dimension, a stronger effect occurs: each rare 



region can independently develop a nonzero steady state 
density. This leads to a smearing of the global transition. 

The smearing mechanism found here relies only on the 
existence of a true phase transition on an isolated rare 
region. It should therefore apply not only to the DP 
universality class, but to an entire family of nonequi- 
librium universality classes for spreading processes and 
reaction-diffusion systems. Note that while the presence 
or absence of smearing is universal in the sense of critical 
phenomena (it depends on symmetries and dimensional- 
ity only), the functional form of the density and other 
observables is not universal, it depends on the details 
of the disorder distribution |3l| . Smearing phenomena 
similar to the one found here can also occur at equilib- 
rium phase transitions. At quantum phase transitions in 
itinerant electron systems, even point-like impurities can 
lead to smearing In contrast, for the classical Ising 
(Heisenberg) universality class, the impurities have to be 
at least 2d (3d) for the transition to be smeared which 
makes the phenomenon less likely to be observed |3^| . 

In conclusion, extended defects destroy the DP transi- 
tion by smearing and lead to a (nonuniversal) exponen- 
tial dependence of the density and other quantities on 
the control parameter. We suggest that this disorder- 
induced smearing may be related to the striking absence 
of DP scaling JTlJ in at least some of the experiments. 
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versity of Missouri Research Board. 
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